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Abstract 

We consider f(1Z) modified gravity theories for describing varying 
speed of light in a spatially flat FRW model, and find some exact solutions. 
Also we examine the dynamics of this model by dynamical system method 
assuming a ACDM background and we find some exact solutions by con- 
sidering the character of critical points of the theory in both formalisms. 
The behaviour of the speed of light is obtained. 

1 Introduction 

The late-time accelerated expansion of the universe is an important challenge 
to the cosmological models. There is an observational evidence based on Type 
la supernovae standard candles [1 and also on standard rulers [H [3] that the 
Universe is in a phase of accelerating expansion now. There are several theoretical 
approaches^ [5] towards the understanding of this accelerating expansion. The 
simplest one is to assume the existence of a positive cosmological constant which 
is small enough to be dominating only at recent times. The predicted cosmic 
history (assuming spatial flatness) is then 

(^) 2 = H%[n Qm (i + zf + nor(i + zf + n A ] (i) 

which provides an excellent fit to the observational data [2]. Models with a 
cosmological constant suffer from the coincidence problem which is the need for 
an extreme fine-tuning of the cosmological constant. To address this problem 
two classes of models have been proposed. In the first class one attributes 
the accelerating expansion to a dark energy (usually the energy of a scalar 
field called quintessence) which has repulsive gravitational properties due to its 
negative pressure [6]. The role of dark energy can also be played by Chaplygin 
gas[7], topological defects [5], holographic dark energy[5], etc. 

The second class of models look for an accelerating expansion via modi- 
fication of general relativity on cosmological scales. Examples of these mod- 
els include scalar-tensor theories [TOI ITT] . f(JVj modified gravity theories [T2]. 
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braneworld models [15]. and so on. These models naturally allow [TU1 [Tl] for a 
super accelerating expansion which the effective dark energy equation of state 
parameter w = p/ p crosses the phantom divide line w = — 1. Such a crossing is 
consistent with some current cosmological data|15). 

The advantage of f(1Z) theories of gravity is that no extra degree of freedom 
is introduced and the accelerating expansion is produced by the Ricci scalar 
(dark gravity) whose physical origin is well understood. This is in contrast 
to other models where the origin and physical meaning of the extra degree of 
freedom is unclear. It has been shown that for appropriate forms of f(1Z) the 
action can produce accelerating expansion at late times in accordance with Snla 
data [TB] . On the other hand, the main drawback of such theories is that they 
are seriously constrained by local gravity experiments [171 [T8l [T9] . 

It can be shown[17] that f(lZ) models are dynamically equivalent to scalar- 
tensor theories with vanishing Brans-Dickc parameter (oj = 0) and a special type 
of potential. This implies that in principle the reconstruction of f(7Z) from a 
particular cosmic history H(z) can be performed in a similar way as in the case 
of the scalar-tensor theories [115] [T3]. However, the vanishing of the Brans-Dicke 
parameter requires some modifications of the reconstruction methods especially 
when the reconstruction extends through the whole cosmic history through the 
radiation and matter eras. The dynamical systems approach followed in the 
present study illustrates these modifications. 

The construction of cosmological models incorporating late accelerating ex- 
pansion based on f(TZ) theories has been an issue of interesting debate recently. 
This debate started from Ref.[2~Q] which demonstrates that f(TZ) theories behav- 
ing as a power of 1Z at large or small 1Z are not cosmologically viable because 
they have the wrong expansion rate in the matter dominated era. This con- 
clusion was challenged in Ref.[5T] claiming that wide classes of f(TZ) gravity 
models including matter and accelerating phases can be phenomenologically re- 
constructed by means of observational data. The debate continued with the 
recent work [55] where a detailed and general dynamical analysis of the cosmo- 
logical evolution of f(lZ) theories was performed. It was shown that even though 
most functional forms of f(1Z) are not cosmologically viable due to the absence 
of the conventional matter era required by data, there are special forms of f(TZ) 
that can be viable with appropriate initial conditions. 

There are some ideas suggesting that the constants of nature, such as the 
speed of light, should be space-time dependent [23] (24] [25] [26]. Theories with 
varying speed of light (VSL) have been firstly proposed by Moffat, Albrecht, 
Magueijo and Barrow [25] [37] as an alternative approach to the inflation mech- 
anism for solving some problems of Big-Bang cosmological models 26, 28]. In 
their formulation the Lorentz invariance is broken and there is a preferred frame, 
in which the speed of light depends only on time. In this frame there exists a 
pre-set function [57] [5U] representing the speed of light and enters in the Fried- 
man equations as an input. 

It is a well-known fact that it is possible to have a varying speed of light 
theory and preserving the general covariance and local Lorentz invariance [3 Oj. 
The price that have to be paid for this, is to introduce a time-like coordinate 
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a; which is not necessarily equal to ct. In terms of x° and x, one has local 
Lorentz invariance and general covariance. The physical time t, can only be 
defined when dx° /c is integrable. 

The most general scalar-tensor action of gravity which allows for a dynamical 
speed of light is illustrated in [31] . This action is previously analysed by many 
authors. Demianski et al. [35] present a class of cosmological models derived 
from Noether symmetry requirement. These models describe accelerating evo- 
lution of an FRW universe filled with matter and exhibit power law dependence 
of the coupling factor and the potential to the scalar field. There is also some 
tracking solutions of this model, in which the time evolution of the scalar field 
tracks the expansion rate of the universe. 

Here we shall investigate the exact cosmological solutions with varying speed 
of light in the framework of f(7V) modified gravity theories. In the following 
sections we shall find some exact cosmological solutions for the spatially flat 
universe. In section 3 we shall examine the dynamics of this theory by dynamical 
system method assuming a ACDM background. Considering the character of the 
critical points of the theory we find some exact cosmological solutions. In section 
4 we shall use the solution of the above system to reconstruct the cosmological 
evolution and functional form of the function f(R) and the speed of light. 



2 The Model 

The action which we use here is similar to the Jordan- VSL action in [31] , except 
that we have changed the gravitational part to f(TV)- 

S = J d*x^(hWf(K) - 2[/(V0 - Z{f)g» v d^d v f) + S m [<f>i;g^] 

(2) 

Here h(ip) — (^-) 4 and U(ip) are arbitrary regular functions of the scalar field ip 
(the field that generates varying speed of light), representing the coupling of the 
scalar field ip with geometry and it's potential energy density respectively. Co 
is the constant velocity of light and hereafter we shall put 8irG = Cq = 1. The 
first part of the above action functional is the gravitational part, including f(7Z) 
and a dynamical term for the velocity of light with arbitrary coupling function 
Z(ip). The latter is the action of the matter fields, (pi, and we assumed that it 
does not involve the scalar field ip, so that the matter is minimally coupled to 
gravity. As emphasized in the introduction, here it is assumed that there is a 
time-like coordinate a; and since dx° /c is not necessarily integrable, it is not 
always possible to define time. It has to be noted that Z{ip) can always be set 
equal to unity by a redefinition of the field ip. Finally it has to be noted that 
here we adopt a metric approach so that the metric, the scalar field ip, and the 
matter fields <pi are dynamical variables. Varying the action with respect to 
metric and ip field gives respectively: 

Kfn^ - lf 9tw ) - V„V„(fc/') + g^V a V a {hf) = 
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Tn„ + d^ipd^ip - ^g^ v (d a tp) 2 - g^U (3) 

and 

„ , dU 1 „dh ... 

V ^ = # - 2 7 # (4) 

where /' is the derivative of / with respect to 1Z. 

The weak equivalence principle holds because the matter fields are minimally 
coupled to the metric. This means that we have V^X 1 ^ 1 = where the energy 
momentum tensor of matter is defined as usual; namely T^ v = -j=^jjn*- In a 
cosmological context, applying the above field equations to the spatially flat 
FRW universe in which the metric has the following form: 

ds 2 = -(dx ) 2 + a(x ) 2 (dr 2 + rW) (5) 

and assuming the matter field as a perfect fluid, we have: 

ihf'H 2 = p + i^ 2 + ±Khf - X -hf + U- 3H(hf>) (6) 

-2hfH = p + p + i, 2 + (hf')-H{hf>) (7) 
1 dh , dU 

^ + 3H ^2dt f -dj (8) 

p + 3H(p + p) =0 (9) 

where a dot over any quantity denotes derivative with respect to the time-like 
coordinate x°. These are c— variable Friedman equations, the field equation of ip 
and the conservation law respectively. H(x°) — is the Hubble parameter, p 
and p are the energy and pressure density of a perfect fluid considered as matter 
field. These equations form a coupled set of non-linear differential equations for 
H(x°) and ip(x ). The time-like coordinate x° is related to cosmic time by the 
relation: 

dt = — = h'Ux (10) 
c 

In the cosmological application dx° /c is integrable and gives the physical time. 
Therefore, the physical Hubble parameter H p (t) = can be evaluated as 
H p (t) = H(x°)*£. Substituting by in eq.@ gives: % + 

3H p (p + p) = This shows that in this model the conservation equation (J9j) is 
valid even in terms of the cosmic time. Since in this model 7Z = 6(H + 2H 2 ), 
we can rewrite the eqs.©, (J7J in these forms: 

ihf'H 2 = p+ i^ 2 + 3(£T + 2H 2 )hf - hif+ 

U — 3H[hf' + 6hf"(H + AHH)] (II) 
-2hf'H = p + p + i) 2 + hf + l2hf"{H + AHH) + 36hf"'(H + AHH) 2 + 



A 



6hf"[H + 4HH] - H[hf + Qhf"{H + AHH)] (12) 



The cosmological solutions of the greatest interest are those for which the 
time evolution of the Hubble parameter is proportional to the inverse of the cos- 
mic time (corresponding to a power law expansion) or a constant (corresponding 
to de-Sitter expansion) . We can thus distinguish two c— dominated uni- 

verse (S m = U = 0), and a (c — A)— dominated universe for which S m = but 
U = Kg is not zero. This corresponds to adding a constant (— 2A) to f(lZ) 
which is equivalent to a cosmological constant. 

In order to have explicit solutions, we have to choose the form of f(TZ). In 
[3"3l the VSL Friedman equations for f(TZ) = 7Z is investigated. But for an 
f(TZ) model with constant speed of light, according to the stability conditions 
for de-Sitter space, we have to choose this form of f(ll): 

f(K)=1l-^+b1l 2 (13) 

Then the condition for the existence of a de-Sitter solution is TZq = %/3/i 2 , while 
the stability condition is satisfied if b > 3 J^ 2 ■ For more details see [34] . Here 
in this section we shall use the above form of / for our VSL f(H) model. 

2.0.1 c dominated universe 

Putting S m = U = and f(K) = K - £ + bV? in the equations ([IT ]) . fTS ]) . 
we get two independent equations. Assuming a power-law dependence for the 
coupling coefficient h(if>), these equations have the following solution: 

H ~ constant, ip ~ e ax ° , h ~ ip 2 (14) 

This is a special choice which is used by many authors [30 , 35 , 36] . The coupling 
function h(ip) ~ ip 2 is a particular case emerged by requiring the existence of 
Ndether symmetry |32] and a is related to the other constants like fi and b. 
The cosmic time and the physical Hubble parameter are: 

a ~ f (15) 

The conditions which one should impose on VSL models are usually inspired by 
the cosmological puzzles. In order to solve the horizon problem of the standard 
cosmology, one should set | — - > for the early universe(see [23] ) and also 
one has a > 0. So requesting an expanding universe together with the horizon 
criteria, one gets the following constraint; v > 0. 

2.0.2 (c — A)— dominated universe 

As mentioned before, this era corresponds to a matter free universe but the 
potential is non-zero. Assuming this form U = Kg in which A is a constant and 
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demanding power-law expansion for cosmic scale factor, one can find that the 
solution is: 

H~ constant, V ~ e ax ° , h^^ 2 . (16) 
and the cosmic time is: 

t-e^ QI °. (17) 



3 Dynamics of VSL f(TZ) Cosmology 

Another way to find out some exact solutions for cosmological models is the 
dynamical system method |37| . In this method by choosing some appropriate 
variables, one can convert the field equations of the desired theory to a set of 
autonomous differential equations. Then the critical points of the autonomous 
system describe interesting exact solutions. Also one can use this method to 
check the stability of the solutions. Dynamics of a scalar tensor theory in 
the Jordan frame using metric approach has been considered in |331 38 . Also 
dynamics of f(lZ) cosmology (not VSL f(TZ) cosmology) has been considered 
in [39]. Here we consider a class of VSL theories described by action ([2]). We 
should notice that the volume element is defined as dx°d 3 x which is different 
from the canonical volume element, and in the field equations H is not the 
physical Hubble parameter and derivatives are with respect to a; coordinate. 

In section 2, we chose h(x°) such that the corresponding solution for cosmic 
scale factor was physically interesting. Here we impose a general form for phys- 
ical Hubble parameter which is related to ACDM cosmology, given by eq.([T]). 
Also we assume U(ijj) ~ h(ip) n where n is a constant [30]. The form of f(lZ) is 
not fixed here. 

Let's rewrite the VSL Friedman equations ([5]), © as: 

3hFH 2 = p m + p rad + \^\? + -HhF - -hf + U- 3HhF - 3HhF (18) 

4 ... ..... 

- 2hFH = p m + -p rad + ip +hF + 2hF + hF - HhF - HhF (19) 

where F := 4L and p m , p ra d represent the matter and radiation energy densities 
which are conserved according to 

p m + 3Hp m = 0, p r ad + p r ad = 0. (20) 

In order to study the cosmological dynamics implied by eas.(|lgp. (IT9| we express 
them as an autonomous system of first order differential equations. To achieve 
this, we first write (|18|) in the dimensionless form as 



, Pm Prad tf 2 K / U h! F' 

3hFH 2 ZhFH 2 6hF 6H 2 6FH 2 3hFH 2 h F ' 

where ' = ^4— = ™ = We now define the dimensionless variables x\,...,X7 

as 

-F' -h' U -f 

x l '■= —^, x 2 ■= —^,X 3 := — ^ TT r, ,X4 ■ 



F h ' 3hFH 2 ' ' 6FH 2 : 
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X5 



n 



H' „ 
— + 2,x 6 



, x 7 := 



Prad 

3hFH 2 



^rad 



6H 2 H ' ' u ' 6/iF' 
As one can see x 7 is in fact fl ra d and cdi + X2 + X3 + 2:4 + £5 + xe 
associated with the curvature dark energy (dark gravity). Defining £l„, 
we can write eq. (l2"T1) as 



(22) 
Ude is 

- Pm 

" 3hFH 2 

Q m — 1 — X\ — X2 — X3 — X4 — X5 — Xq — X 7 . (23) 

Using the defined dimensionless variables, we can express eg. (If 91) as 

x' 1 +x' 2 = — 1 — 3^3 — 3x4 — X5+ 3xq + X7 + x1 + x 2 — X2X5 + 2xiX2 — X1X5. (24) 
Also, differentiating X3, ...x 7 with respect to N we have 

x 3 = ^[(l - n)x 2 +x\ — 2x 5 + 4] 

X1X5 ,„ .s 

2:4(2x5 — xi — 4) 



m 



X1X5 



in 



2x 5 (x 5 - 2) 



x 6 (-2 
x 



where 



2x 5 + xi + x 2 ) + x 2 (x 4 + nx 3 ) 
7 - x 7 (x! + x 2 - 2x 5 ) 



(25) 
(26) 

(27) 

(28) 
(29) 

(30) 



/' f.TZ 

and ^ implies derivative with respect to 1Z. 

The autonomous dynamical system ([2"3j). PS). pj]). (|2"7) . (l2"5)l and K2SJ) is 
the general dynamical system that describes the cosmological dynamics of VSL 
f(TZ) theories. Instead of investigating the above autonomous system for various 
different behaviours of m(f(7i)) we eliminate m from the system by assuming 
a particular form for H p {N) consistent with cosmological observations. Once 
Xf,{N) is known we can solve (|2T1) for and substituting in (f2l)|). we find 



= — Xg — 2X5(X5 — 2) — X4(2X5 — Xi — 4) 



(31) 



which along with (|2"4"|) . (|2"5|) . (|2"8)l and (|2T))) describes a new dynamical system 
which is independent of m. On the other hand one can easily verify that = 

\x2-, so by substituting this relation in the above equations, we have: 



-3— 3x 3 — 3x 4 — — -(l+xi+x 2 ) — ^x 2 +3x 6 +x 7 +x^- 



-Xi 



7 



7X2 



x 3 



1 H' 

{--n)x 2 +xi-2j^ 



H' 1 



x' = -2 + ±x 2 + 2 -£ + ± : 



i? 



^ 1 



r x 2 



X4 



H' 
2-?- 

H„ 



2 X2 



Xl 



(32) 
(33) 

(34) 
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(H f 1 \ 

-6 - 2j^- + x\ + -^ 2 J + x 2 (^4 + (35) 

(1 i?' \ 

an + -x 2 ~ 2 Jf-4j (36) 

and to = H r Xl — . 

The results of our analysis do not rely on the use of any particular form of 
0:5 (JV) (i.e. H p (z)). They only require that the universe goes through the radia- 
tion era, matter era and acceleration era. For the sake of definitcness however, 
we will assume a specific form for H(z) corresponding to a ACDM cosmology 
([T]l which in terms of N, takes the form 

H p (N) 2 = H$[n 0m e- 3N + n 0rad e- 4N + Q A ] (37) 

where N = In a = — ln(l + z) and Q\ = 1 — Qo m — Qorad- It is straightforward 
to study the dynamics of the system (f32|) . (|33]l. (|34|) . (|35|) and ([36]) by setting 
= to find the critical points and their stability in each one of the three eras. 
By setting x[ = 0, we have 

H 1 1 

£5 = tt + 7 a; 2 + 2 (38) 
Tip 4 

These equations describe the cosmological dynamics of the VSL f(lZ) theory. 

In general, we should add a further parameter which could be related to 
and h^, beside the above n. In fact if we do not try to reconstruct 
the function h(ip), such a function can be fixed priori and the corresponding 
parameter would be, for example, In such a case H{N) would not be 

fixed as in our reconstruction approach but would have to be determined by the 
autonomous system. 

The critical points are shown in Tables (1), (2) and (3). The stability anal- 
ysis of these Tables assumes that X5 = const and therefore it is not identical to 
the full stability analysis where 25 would be allowed to vary. The usual stability 
analysis of the cosmological dynamical systems assumes a particular cosmologi- 
cal model (e.g. a form of f(1Z) or to or n) and in the context of this physical law, 
the stability of cosmic histories H(N) is investigated. In this context clearly 
a stable cosmic history is the one preferred by model. In the reconstruction 
approach, however, the stability analysis has a very different meaning. Here 
we do not fix the model h(i(>) and f(TZ) (physical law). Here we fix the cosmic 
history and allow the physical law f(lZ) and h(ip) to vary in order to predict 
the required cosmic history. Thus our stability analysis concerns the physical 
law f(7Z) and h(tp) and not the particular cosmic history. The physically in- 
teresting quantities are the values of the critical points we find in each era in 
the context of the ACDM cosmic history. These tell us the possible physical 
laws f(TZ) and h{ip) that can reproduce a ACDM cosmic history. Some criti- 
cal points in each era are not stable. This however does not imply that these 
points are not cosmologically relevant. These instabilities are not instabilities 
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of the trajectory H(N) (which we keep fixed) but of the forms of f(Ti) which 
is allowed to vary. Thus they are not so relevant physically since in a physical 
context f(TZ) is assumed to be fixed a priori. Calculation of eigenvalues are too 
long, but in what follows we don't need them. So we will not write them, and 
thus we will not present them here. 



Xi 



X2 



%3 



x 5 



x e 



X 7 







-12 



4m 
1— m 



-8m 
L — m 



2m 



I— m 



(l-m)(l-m a ) 



(l-m)(l-m 2 ) 



8m(n-l)~ 
-l+2n+m 



8(m+l) 
-l+2«+m 



■I ? 



-l+2n+m 



2(m+l) 
-l+2ra+m 



m(4-i^) 



4(l+m) - 



Table 1: The critical points of the system in the radiation era. 

In Table (1) the critical points of the system for radiation era is given, in 
which 

t _ -3(1 - to) + 2n(l - 2m) 

X3 ~ n(-l + 2n + to) 2 (12n(l + to) + 9(1 - to) - 6n(l - 2to)) 

((—1 + 2ra + to) (n 2 + (n — 3)(n - 2) - TOn(53 - 16to 2 + 8m(l + An))) 

-16n(l + to) 2 (to - 2) (to - 3/2)) + — 1 

v ' K A 1 " n(-l + 2n + m) 



and 



* = 4(1 + to) 

X(i ~ (-1 + 2n + to) 2 (12n(l + to) + 9(1 - to) - 6n(l - 2m)) X 

((-1 + 2n + to) (n 2 + (n - 3)(n - 2) - mn(53 - 16m 2 + 8m(l + 4n))) - 



and 



16n(l + to) 2 (to - 2)(m - 3/2)) 



(1 - 16m 2 ) -2 - 27to - 12to 2 + 16to 3 „ 1 2 7 



+ 



3 12(1 + to) ' z 12 

and x\ satifies a cubic equation whose solutions are 

, _ 4(-1 + 2to) 

TO - 1 

„ _ -12to 2 + 16to 3 - 2 - 27to + V64m 4 - 39to 2 - 80to 3 + 52 + 148m 
2 2to 3 — 5to 2 — to + 6 
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Xi 


2m 
1—m 





3m 
1—m 


6m(- 1+nJ 
-l+2ra+m 


m(3 — \x 2 ) 


X2 


2(l-3m) 


6 


— 6m 
1—m 


6(l+m) 


x\ 


1—m 


-l+2ra+m 


X3 











t 





x 4 


l-2m 





-l+4m 


n+2m+l 


3 ,1T 1_T^ 
2 ' 2 2 8 2 


1-m 2 


2(l-m 2 ) 


2n+2nm+m 2 — 1 


(l+m)(-3+i^) 


x 5 


I — 'Am 


2 


1— 4m 




4 X 2 ^ 2 


1—m 


2(l-m) 


-l+2ri+m 


x 6 


6m 2 — 5m+l 


-7 


(1— Am)m 


4 


r tt 


(1—m) (1—m' 2 ) 


(l-m) 2 (l+m) 


x 7 


mibm 2, — 10m+3 
(1—m) (1—m 2 ) 















Table 2: The critical points of the system in the matter era. 



„ _ -12m 2 + 16m 3 - 2 - 27m - V64m 4 - 39m 2 - 80m 3 + 52 + 148m 
2 2m 3 — 5m 2 — m + 6 

In Table (2) the critical points of the system for the matter era is given, in 
which 

t 1 

x\ = x 

3 1 + ra 

6n 2 m 3 - n 2 m 2 - 7mn 2 + 2n 2 - 5nm 3 + 7nm 2 + 8nm - 8n + 3m + 8 - 4m 2 - m 3 
4n 2 + 4nm — 4ra + m 2 — 2m + 1 

and 

xt ■ 



6n 2 m 2 + hmn 2 — 2n 2 + nm 2 + 7nm + 8n + 2m + 1 



4n 2 + 4nm — An + m 2 — 2m + 1 
and 

tt = 2 l(l+m)(-3+^x*) 

X " ~ (l-m)(3-|4) 
and x\ satisfies a cubic equation, whose solutions are 

I 6m 



"2 -, 
m — 1 



2 (6m 3 - 8m 2 - 5m + 6 + V49m 4 - 76m 3 - 44m 2 + 72m) 
—5m 2 + 2m 3 + 6 — m 



t 2 (6m 3 - 8m 2 - 5m + 6 - \/49m 4 - 76m 3 - 44m 2 + 72m) 

2 —5m 2 + 2m 3 + 6 — m 

Finally the critical points of A— dominated era is given in Table 3, where 
x\ satisfies a cubic equation, whose solutions are 

n _ - 6 

,h 9 — 



b 2 — 7 
m — 1 



k _ 2 (5m 2 - 7m - 9 - V49m 4 - 154m 3 + 7m 2 + 210m + 9) 
2 —5m 2 + 2m 3 — m + 6 
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Xi 





—4m 
1 — m 


-1 q 





X2 





8 

1 — m 







x 3 











l—m 
l+m 




-1 


-4+2m 
1 — m 2 


(l+m) ( 2 + 4 X 2) 


-2 
l+m 


x 5 


2 


4-2m 
1 — m 


42-2 + 2 


2 


x 6 





4(-4+2m) 







(l-m)(l-m 2 ) 


(l+m)(-6+^(l-m)xJ) 


x 7 





Q(l-m)+rrv £ (13-5m) 
(l-m)(l-m 2 ) 









Table 3: The critical points of the system in the de-Sitter era. 



^ 2 (5m 2 - 7m - 9 + V49m 4 - 154m 3 + 7m 2 + 210m + 9) 
—5m 2 + 2m 3 — m + 6 

4 Reconstruction of f(H), h(i/j), U(iJj) 

We can now reconstruct the form of the function f{R) and h(ip) and U(ip) that 
correspond to each one of the critical points of the system shown in tables (1), 
(2) and (3). These reconstructions are effectively an approximation of these 
functions in the neighbourhood of each critical point. 

Consider a critical point of the form {x\, X2, X3, £4, X5, xq, £7). Now we want 
to reconstruct the form of the potential U(ip) and h(ip). Using (|22"]) . we find 
that: 

F(N) = F a e- SlN , h(N) = h e~ S2N (39) 

where Fq and ho are the present value of F and h. Assuming that at the present 
time the velocity of light is equal to Co we have to set ho = 1. Using eq. (|22l) 
we find 

1 

MN) = -2{6F Q )i X s e - (gl+S2) ^ + C (40) 
Xi + x 2 

where C is a constant which can be put equal to zero by a shift in ^-field. 
Equations (f3T)f and P0"|) allow us to eliminate N in favour of tp 



h(ip) = £ip (*i+»3) (41) 



where g = (Sl t g2 



It is interesting to note that for all solutions obtained here from the analysis of 
the critical points and for any epoch, the h function has a unique form. Since 
c/co = h 1 ^ we have: 

- = a^ 4 (42) 



Co 

Depending on the value of x-i given in tables (1), (2) and (3), this can lead to a 
constant, decreasing or increasing speed of light with respect to the scale factor. 
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Using the horizon criteria one can obtain some restrictions on xi and thus on m. 
For radiation era this condition leads to X2 > 4. For matter era we have 22 > 2. 
And finally for de-Sitter era it is X2 > —4. This can be used as a selection rule 
for model parameters m and n. 

In a similar way we can reconstruct U(ip). From eq. (|22p . we have 

U(N) = 3x 3 h{N)F(N)H(x ) 2 = 3x 3 h(N)i F{N)H p (N) 2 (43) 

Using now the input form of H P (N) feq. (|3"T)) ), we find the dominant term of 
H P (N) in each era. Using eq.fljD} and (|%3"jl we have: 

U{iji) = X^^rk (44) 

where 

A = 3x 3 Fo^lor£, * 2 i Rod. era 

= 3x 3 F ^om? 13 , Mat. era 
= 3x 3 F a ^ + ^ (1 - fl 0r - Oom), d.S. era 

We can now reconstruct the form of the function f(TZ) corresponding to 
each one of the critical points of the system. This reconstruction is effectively 
an approximation of f(TZ) in the neighbourhood of each critical point. It is 
particularly useful because most of the dynamical evolution takes place close to 
fixed points. Using the relation 

n = 6{2H 2 + H'H) = 6 (2HI fe-i + H' p H p h-^ - -H^'h'^ (45) 

In terms of H p in each era, we obtain 

U(N) = ^ 2 ^0re (i *" 4)Ar , Rad. era 



('+!) 



^n„,e ( ^" 3)JV , Mat. era 



which leads to 



6(2 + ^)(1 - Q 0r - n 0m )e^ X2N , d.S. era (46) 



= F °l.6(2 + f)(i-n„ r -t!„j.l " ' d - s - era (47) 
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and by integration we get 



/(ft) ~ ft -*+*•» 



+1 



Rad. era 



ft 




+1 



Mai. era 




d.iS. era 



(48) 



Although we can reconstruct /(ft) for any critical point and for each one of 
the three epochs in terms of the values of x\ and X2, it has to be noted that all 
of these can be written as 



as it is clear from the definition of m in equation (f30|) . 

5 Conclusion 

Here we have investigated analytically the behaviour of VSL /(ft) gravity. We 
saw that it is possible to produce the background expansion history H (z) in- 
dicated by observations. Choosing the form /(ft) ~ ft 1+m we get c ~ a~ X2 ^ 4 . 
The horizon criteria puts some limitation on x-i and thus on m. It has to be 
noted that the form of the potential of the VSL field has to be fixed according 
to equation (|44[) . 
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